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Abstract. A fourth order nonlinear evolution equation is derived from a mi- 
croscopic model for surface diffusion, namely, the continuum solid-on-solid 
model. We use the method developed by Varadhan for the computation of the 
hydrodynamic scaling limit of nongradient models. What distinguishes our 
model from other models discussed so far is the presence of two conservation 
laws for the dynamics in a nonperiodic box and the complex dynamics that is 
not nearest-neighbor. Along the way, a few steps has to be adapted to our new 
context. As a byproduct of our main result we also derive the hydrodynamic 
scaling limit of a perturbation of continuum solid-on-solid model, a model that 
incorporates both surface diffusion and surface electromigration. 
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1. Introduction 

A process of great technological importance, molecular beam epitaxy (MBE), 
is used to manufacture computer chips and semiconductor devices, see Barabasi 
and Stanley In general the chip is constructed by spraying a beam of atoms 
on a flat surface. There are three phenomena that take place in the construction 
process namely deposition, diffusion and desorption. The atoms arrive or deposit 
on the surface and do not stick on the first contact point, but diffuse or walk on 
the surface. When an atom reaches the edge of another wandering atom, the two 
atoms meet or glue together, forming islands. Sometimes an atom may jump out of 
the surface. Smaller islands may develop into larger islands affecting the roughness 
of the surface on the macroscopic scale. A rough surface does not have very good 
contact properties and engineers would like to understand the basics mechanisms 
affecting the morphology in general. 

How deposition and desorption affect the morphology of the surface, it is quite 
well understood. It is of great importance to know how the profile of the surface 
evolves on the macroscopic scale if the atoms that make up the surface diffuse. We 
assume that no atoms arrive or leave the surface. 

The present paper discusses a model for surface diffusion the continuum solid-on- 
solid model, known also as the fourth-order Ginzburg-Landau model. The system 
has a complex interaction that is not a nearest-neighbor interaction and has two 
conserved quantities in a nonperiodic box, namely the total slope of the surface 
and the linear mean of the surface slope. Under the assumption that the molecules 
of the surface follow the dynamics of the continuum solid-on-solid model, we will 
prove that the dynamics of the surface slope profile on the macroscale is a fourth- 
order nonlinear equation. Since the model is nongradient, the derivation of the 
limit is not trivial and we shall use the method developed in Varadhan |17j and 
further extended in Quastel Varadhan and Yau 18 . Even though we have 
an interaction that is not of nearest-neighbor type, the microscopic current still 
decomposes as the Laplacian of the slope field and the fluctuations. Usually in 
nearest-neighbor models after subtracting the fluctuations from the current we end 
up with the gradient of some local function. We use the result that for continuum 
solid-on-solid model the space of exact functions has codimension one inside the 
space of closed functions, see Savu [T^ . 

We also include a discussion of a perturbation of continuum solid-on-solid model, 
where the evolution of the surface is driven by both diffusion and electric field. The 
electric field will add one extra second-order term to the final nonlinear equation. 

The continuum solid-on-solid model is an approximation of the discrete solid-on- 
solid model and hence is not considered a truly microscopic model. Unfortunately 
at the time of writing of this paper, we don't have the required techniques to solve 
the discrete solid-on-solid model. However, partial rigorous results and numerical 
analysis are available for the discrete case, see Krug, Doobs and Majaniemi 

Finally the paper is organized as follows: section 2 contains the description 
of continuum solid-on-solid model and of the model for surface electromigration, 
the statement of the main results and a note on similar models considered in the 
literature so far; section 3 outlines the proof of the main result, the computation of 
the scaling limit of continuum solid-on-solid model; section 4 shows how the final 
nonlinear equation is identified; in section 5 we calculate the asymptotics of central 
limit variance, the main ingredient used in section 4 to prove that the microscopic 
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current can be replaced by a multiple of the field Laplacian; in section 6 we calculate 
the scaling limit of the modified model to incorporate surface electromigration. 

2. The models 

The continuum solid-on-solid model, is a 1-dimensional lattice system with con- 
tinuous order parameter, used to study the evolution of an interface. The model 
describes the movement of the interface at the mesoscopic level and hence is not 
considered a truly microscopic model that captures all the aspects of the particle 
motion, but it has the advantage of being more suitable for computations. 

The height model. Let T be the torus represented as the interval [0, 1] with 

and 1 identified. For each positive integer N, there arc scaled periodic lattice 
points located at sites i/N inT, i = 1, . . . ,N. We shall denote by hi{t) the height of 
the surface at the site i/N, at time t. Also because of the periodicity of the lattice 
points, hN+i{t) = hi{t). The energy function H[^{h) of a height configuration h is 
chosen to be invariant under the global translation hi{t) i — > hi{t) + c and has the 
form 

N 

HN{h)=Y,V{hi+i-hi). 

i=l 

In this paper we shall assume that the potential is quadratic V{x) = x^. We shall 
require the evolution in time of the surface to preserve the sum of the heights, to 
have as invariant measure the infinite mass measure e^^'^'^^^dh, and to be reversible. 
All these properties are satisfied by the solution of the stochastic differential system: 

dhi{t) = — —{wi - Wi-i)dt + N^iy/aTidBi - ^amdSj.i), 1 < i < AT. 

Above, ai{h) = a{hi—hi-i,hi+\—hi, hi+2—hi+i) where a is a function with bounded 
continuous first derivatives that satisfies < 1/a* < a(.T_i, .tq, < a* < oo. The 
n copies of the Brownian motion Bi, i = 1, . . . ,n are independent. We define the 
instantaneous current 'Wi{h) of particles over the bond 

Wi{h) = {d-ia-2doa + dia){hi-hi-i,hi+i-hi,hi+2-hi+i) + 

(2.1) +ai{h){V'{hi - hi_i) - 2V'{hi+^ - hi) + V'{hi+2 - hi+i)). 

Because the height process preserves the sum of the heights of the surface, this 
dynamics models surface diffusion. 

The slope process. A crucial property of the dynamics of the heights is the 
gauge property, namely the dynamics invariance under the action of the group G 
of translation in the (1, . . . , 1) direction, 

G = {T : ^ I T(a;i,...,a;jv) = {xx + c, . . . ,xn + c), c e M}. 

Hence there exists an induced dynamics on the quotient space R-'^/G of equivalence 
classes. A representative of an equivalence class is the slope configuration Xi{t) = 
h^+i{t)-h^{t), l<i< N-limdxN{t) = hi{t)-hN{t). Note that J2^^l^^it) = 0- 
As a function of the slope configuration of the surface, the energy becomes Hj^{x) = 

EliVix,). 

In the sequel, we shall study the slope process rather than the height process. 
The slope process is reversible and has as equilibrium distribution, the product 
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probability measure dv'j^ — e ^"^^^^/Z^rfx. Below we write down the stochastic 
differential system, the generator, and the Dirichlct form of the slope process: 
(2.2) 

TV* 

dxi{t) = — (wi+i - 2mj + Wi^i)dt + N'^{,Ja~[dBi^i - 2^/a7dB,; + ^a^+idBi+i), 



(2.3) N^LN{f) - — ^a.(a:)(a.+i- 29.+ 5.-1)2/ + w;,(a:)(a.+i- 29,+ 5.-i)/, 

AJ-i f JL p-Hn(x) 

(2.4) N^DMif) = — V a,{x){d,+^f - 2d,f - d,-if)\x) dx. 

The factor in the generator with the lattice spacing of 1 jN represents the scaling 
of space and time. This scaling is needed to observe a nontrivial motion in the limit. 

The diffusion (|2.2() is driven in the direction of the linear vector fields Xi — di+i — 
2di + 9i_i, I < i < N, therefore is not ergodic in the whole space M^. It becomes 
ergodic when restricted to the hyperplane xi + ■ ■ ■ + xn = Nx of average slope 
X. The unique equilibrium probability measure of the dynamics restricted to the 
hyperplane is the conditional probability e^^'^''^'^ / dx given xi + • • ■ + xi^ = Nx. 

Instantaneously, the slope profile of the surface decreases at some site i twice as 
much as increases at the adjacent sites i — 1, i + 1. We see that any update of the 
slope configuration affects the slopes at three sites. This type of interaction, known 
as the three site interaction, is quite complex and has been very rarely studied, so 
far. 

We should note an integration by parts property of the current Wi, for all bounded 
and smooth functions /: 

(2.5) E'^[w, ■ /] = E'^[a, ■ (a,+i - 2d, + a,_i)/]. 

The expected value above is with respect to the equilibrium measure dv'^. 

The dynamics in a nonperiodic box. Although we are concerned with the 
study of the slope process defined on a periodic space, we make use of a similar 
slope process evolving in a box with nonperiodic boundary. Below we describe this 
new dynamics and its main properties. Suppose A is a box with finitely many sites 
of the lattice Z. The infinitesimal generators, 

(2.6) La(/) = i a,ix){d,+i-2d,+d,-iff+w,{x){d,+i^2d,+d,-i)f, 

ieA,z+l,i-16A 

respectively, 

(2.7) L^(/) - 2d, + a,_i)2/ + w,{x){d,+i - 2d, + 9._i)/, 

produce two diffusion processes. 

In the case of the first dynamics, generated by La, there is no transport of par- 
ticles over the boundary of the box, so the box does not have periodic boundary. 
If the box A = [—1,1] D Z or A = [i — l,i + I] D Z we use the shorter notation L;, 
respectively, Li^i for the generator La- The second dynamics generated by Loo is a 
dynamics on the infinite lattice Z. The dynamics Li^i preserves the average slope, 
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yli = ' 2^;+-i"^'^' ; the linear mean of the slopes, yfi = ^ ; inside 

the box [i — l,i + ^D'Z. The second conserved quantity, yf[, should be understood 
as a boundary condition that is preserved in time. The linear mean slope, 
is conserved in time because we have a model for surface diffusion and the total 
height of the surface does not changes as time passes by. We will use two different 
notation Xij and y^i for the mean slope of the field in the box centered at i, of 
size I. Also yij stands for the vector of the conserved quantities, {yli,yfi). As a 
convention we will drop the subscript i, meaning the center of the box, if the box 
is centered at the origin. 

Equilibrium measures. We proceed to describe next the equilibrium measures 
of the dynamics that we have introduced. For the dynamics the grand canonical 

measure is the product probability measure v^'^^ ^ — - — ^^^^ ' dxj, whereas 
the canonical measure is the conditional probability measure i^y n = i'^^^ ;( \yi,i) 
given the level set 

^"^"^^ I 2]T1 liHT) 

The canonical measure is the unique stationary probability measure for the re- 
stricted dynamics on this set. The Dirichlet forms of the operators Li, Lu, with 
respect to the grand canonical measure with a — 0, are denoted by Di{f), re- 
spectively, Di^i{f). The Dirichlet form of the operator Li with respect to the 
canonical measure is denoted by D^n ^{f). For Laa, the product measures j/^^ — 

~ — ^z{a) ' ^'^^ equilibrium measures. We shall use the notation for the 

product probability measure ~ — -^-^dxi. 

A model for surface electromigration. We consider also a perturbation of 
the continuum solid-on-solid model. The new system describes the evolution of a 
one-dimensional surface driven by both surface diffusion and surface electromigra- 
tion. The surface electromigration refers to the motion of atoms on a solid surface 
that is caused by an electric current in the material. The electric field interacts 
with the atoms of the surface as the wind blows the sand particles, and a ripple 
pattern is observed in the long run. Electromigration along interfaces is believed 
to play a crucial role in the failure of metallic circuits, see Schimschak and Krug 
[TS]. for further details. 

We assume the electric field is a continuous function E{t, 6) defined on [0, T] x T. 
As before Xi represents the slope of the surface at the site i/N. The generator of 
the system on a periodic lattice, that incorporates the action of the electric field, 
is 

(2.8) 

N^Ln.eU) - N^L^if) + —Y,E[t,-ja{x,^l,x,,x,+l){^^^l - 2d, + d,+i)f. 

i=i ^ ^ 

Hydrodynamic scaling limit of the models. We shall call ^'^"^ -^i?t 

the law up to time T of the slope process (12. 2f) started in some nonequilibrium 
, equilibrium distribution z^^, respectively. The law up to time T 
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of the slope process (|2.8ll , driven by the electric field, started in the nonequilibrium 
measure i''^'^ shall be called -p- 

Under P^"^ and P^rp the random variable 

(2.9) TTNit) = ^ + ■■■ + XN{t)Sr 

has distributions and rp, respectively. We also refer to the random variable 

(|2.9(l as the empirical distribution. Every realization of this random variable is a 
measure-valued continuous path and Q^^\ is a distribution on the space X — 
U(;>o) ^([0; T],Mi). The space X is endowed with the inductive hmit topology, the 
strongest topology that makes all the inclusions of C([0, r],7Vl;) continuous. The 
space of signed measures A^;, with total variation not exceeding is a metrizable 
space with the weak topology. 

We say that a model has hydrodynamic scaling limit if under certain assumptions 
the sequence of laws of empirical distributions has a limit that is supported on the 
solution of an initial value problem. 

Definition 2.1. A sequence of initial distributions u^"^ onM.^ is said to correspond 
to the macroscopic slope profile tuq G L^{T) if the random variable ttjv converges 
weakly in probability to 5mo{e)de, i-^-j for any continuous function (f> G C(T) and 
any e > 0, 



1 ■ ncq 
iim sup Vj^ 



N , . 
1 \ ^ , / I 



2=1 



(j){e)moi9)d0 



T 



>e'> =0. 



Our first result of the paper says that the continuum solid-on solid model has a 
scaling limit and provide the form of the limiting evolution equation. More precisely, 

Theorem 2.1. Assume that the potential V{x) is equal to x^/2. Let mo e L^(T) 
be a macroscopic slope profile such that Jj.mo{0)dO = 0. Assume that the sequence 
of initial distributions {^''j^'^}n corresponds to the profile TOq and the initial relative 
entropy H is of order 0{N). Then the sequence of probability measures 

!N,TS 



{Q^'t} N>o is tight in X . 



Any possible limit, Qt, of a convergent subsequence of {Q^'^rp} n>o is concen- 
trated on the weak solutions {m{t,6)d9)ffz[Q^T] of the Cauchy problem with periodic 
boundary conditions 

(2.10) dtm = -^d^{a{m)d'^m), m(0, 9) = mo{9), 9 eT. 

The transport coefficient d is a nonrandom continuous function on M., given by the 
following variational formula, 

(2.11) a{a) = inf ' ' 
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a{x-i,xo,xi) ( 1 + - 2di + 5j_i) ^ r-'g^ j 



The infimum on the line above is taken over all local functions g{x^s, ■ ■ ■ , Xs, Vq i) 
of the slope configuration. The shift acts on configurations {T^x)k — Xk+j and 
on local functions, {t^ g){x) = g^r^x). The expectation E^a^ is with respect to the 
grand canonical measure . 
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Note. It can be shown that the transport coefhcient a is a continuous, bounded 
above and below function (see Kipnis and Landim [H]). Also we expect that the 
methods of Landim, OUa and Varadhan |1(J| can show that a is smooth but we will 
not pursue it here. 

Note. By a weak solution of the Cauchy problem (|2.10() we mean a path m ^ X 
such that for each time T > 0, the value of the path m at the time T is a Lebesgue 
absolutely continuous measures on T, satisfying the energy estimate 

(2.12) ( [ {d^m{t, 6) fdOdt <oo. 

Jo Jt 

Moreover for each < T < cx3 and for each test function (f) e C'^''^{[0, T] x T), 

[ m{T,e)(p{T,e)d0 ^ f m{O,e)(p{O,9)d0 ~ ( [ m{.s,e)d,<j){s,9)deds+ 
Jt Jt Jo Jt 

+ 1 [ [ a{m)d^m{s,9)d^(j){s,9)d0ds ^ 0. 
2 Jo Jt 

If we assume that the initial condition toq has the property that J,j,'mo{9)d9 — 
then for each time t > the solution satisfies m{t, 9)d9 = 0. 

Uniqueness of weak solutions of the Cauchy problem H2.10|l . that satisfy the 
energy estimate (|2.12|) has not been proved yet. If the transport coefficient d does 
not depend on the field m, the uniqueness of the Cauchy problem is known and can 
be found in Eidelman book j4|. 

As will be explained later the fluctuation dissipation equation for the continuum 
solid-on-solid model follows from the direct sum decomposition of a Hilbert space 
to be defined next. 

We define the Hilbert space of closed functions, Cx to be the space of those 
^ S L^{dv^) that satisfy in the weak sense the equations Xiir^S) — Xj{T^^) for 
all integers i and j. It is not hard to see that a subspace of Cx is the closed 
linear span in L^{dv^) of functions = T'' 9)1 where 5 is a bounded local 

function with bounded first derivatives. Even though the infinite sum 
does not make sense, the function is well defined because the vector field kills 
all but finitely many terms of the infinite sum. We shall call this space the space 
of exact functions and we shall use the notation £x . The space of exact functions 
has codimension one inside the space of closed functions. In this paper we do not 
include the proof of this result, since it is very technical and is not of probabilistic 
nature, however we include the statement, see Lemma l2. II Lemma l2. II is discussed 
in Savu 

Lemma 2.1. Let 1 denote the constant function 1. The direct sum decomposition 
holds: 

(2.13) Cjf-R10fx. 

The second result of the paper proves that the model for surface electromigration 
(|2.8|) has a scaling limit as well, and calculates the limiting evolution equation. 

Theorem 2.2. Suppose the hypothesis of Theorem \2.1\ are satisfied. Then the 
sequence of probability measures {Qjv b t\n>o is tight in X , and any possible limit 
Qe,t is supported on the weak solutions of the Cauchy problem 

(2.14) dtm^ ~]^dl{d{m){dlm + E)), m{0,9) = raQ{9), 9 eT. 
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The transport coefficient a is given by the same variational formula as in the state- 
ment of Theorem \2.1\ 

Similar models. The continuum solid-on-solid model belongs to a large class 
of Ginzburg-Landau models. The slope model is of nongradient type and has an 
unusual dynamics because two neighboring exchanges occur always simultaneously. 
Nishikawa 11 , and Bertini, OUa and Landim 2 have investigated the hydrody- 
namic scaling limit in higher dimension of the gradient version (i.e., a is a constant 
function) of the slope model. They have found that on the macroscale the interface 
follows a fourth-order nonlinear evolution equation. 

Another similar model, the second-order Ginzburg-Landau model, where the 
sum of the heights is not conserved, was the subject of extensive discussions in 
the literature: the hydrodynamic scaling limit was derived by Fritz jHi and Guo, 
Papanicolau and Varadhan 7 for the gradient version, and by Varadhan for 
the nongradient case, whereas the nonequilibrium fluctuations have been proved by 
Chang and Yau 0. The second-order Ginzburg-Landau model and the continuum 
solid-on-solid model correspond to Glauber, respectively, Kawasaki dynamics in the 
context of interacting particle systems. As expected, a different dynamics at the 
mesoscopic level causes different dynamics at the macroscopic level, a second-order 
parabolic differential equation in the case of second-order Ginzburg-Landau model 
versus a fourth-order parabolic differential equation for the continuum solid-on-solid 
model. 



3. Hydrodynamic scaling limit of continuum solid-on-solid model 

In this section we give a sketch of the main result. Theorem 12.11 We follow a 
standard scheme to derive the hydrodynamic scaling limit of the continuum solid- 
on-solid model. The existence of the limit follows from the tightness of the sequence 
of probability measures {Q^Wn- Let Qt be the limit of some weakly convergent 
subsequence of {Q^%}n- We proceed to characterize the limit Qt, showing that 
it is supported on continuous paths with certain regularity property, known as the 
energy estimate (|2.12ll . The most involved part of the argument is the identification 
of the possible weak limit Qt as some probability measure supported on the weak 
solution of the Cauchy problem (|2.10(l . We shall make the assumption that the 
sequence of initial distributions {v^'^}n corresponds to some macroscopic profile 
Too G L^{T). 

Note on notations. Throughout the paper we shall make use of the shorter no- 
tation limsup^^^j^ Zn^i„ /(-^i, ■ ■ ■ I ^n) for the sequence of limits 
limsup_jj^^j^ . . . Iimsup2^_^j^^ f{zi, dots, Zn). 

Tightness. As a consequence of Prohorov theorem and Arzela-Ascoli theorem, 
the tightness of the sequence follows from the next two Lemmas. 

Lemma 3.1. For any test function (j) G C^(T), any finite time T , and any e > 0, 
(3.1) 

limsup P^";^! sup ^Zl<^( ^ j^^'W-^Zl'^l ^)^»(*) > =0- 
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(3.2) 



r 1 ^ 

limsup P^°^< sup — V 

~+oo,JV-*oo ' L 0<t<T JV ^ 



\x,{t)\ >l\=0. 



Proof. The first convergence H3.1|l follows from Garsia-Rodemich-Rumsey inequal- 
ity. To prove the second convergence (|3.2() we can use estimates on the moment 
generating function of hitting time of the diffusion process H2.2|l . The reader may 
consult Kipnis and Landim 8^ or Guo, Papanicolau and Varadhan or Savu |14| 
for a complete proof. 

It is interesting to note that the stronger superexponential estimates can be 
established for the process in equilibrium, 
(3.3) 



lim sup 



(3.4) 



— loff P ^ 



sup 

0<t,s<T.t-s<S 



1 r 1 ^ 1 

lim sup — log <^ sup —y2\xt{t)\>l\ 

^oo.A'^oo JV 0<t<T ^ J 



w)(/)( — ]du 



> e 



~oo. 



□ 

Energy estimate. Any limiting point Qt of the measure- valued sequence {Q^^jjv 
is supported on paths ^ £ X such that at each time t, fi(t) is a Lebesgue absolutely 
continuous measure on the torus T with density m{t,6). Moreover for each finite 
time T, the density m(t, 6) satisfies the energy estimate 

(3.5) / [ [dlm{t,9)fdedt <oo. 

Jo Jt 

We note that the energy estimate (j3.5(l is equivalent to the inequality 

(3.6) sup / [ {2m{t,e)d^(l)-C(l)^)d0dt < oo, 



sup / 

0eCi-2([O,T]xT) Jo 

where C is a constant not depending on (j). 

The entropy inequality H4.10|l . Feynman-Kac formula H4.12|l . and Lemma l¥!^ can 
be used to derive the estimate 



lim sup sup E'^°'^ 



1 ^ 

N ^ 



2<p[t,-]N^Av^^J-\^,{Ax)r 



(3.7) 



-Vu(Ax,y)02( )df 



< oo. 



Here, the cylinder function Aa; is the discrete Laplacian of the slope field, x_i — 
22:0 + a;i, and the variance Yi^i{Ax, y) is defined later in section 4, see definition ??. 
As will be proved in Lemma l5.2l the variance, Yi^i{Ax, y), has a uniform- in- ?; lower 
bound, therefore the estimate H3.6|) follows from H3.7|l . after integrating by parts. 
Moreover, we can conclude that at any time t the weak second derivative of the 
measure /z(t) is in L'^{T,dO) and hence ^{t) is absolutely continuous with respect 
to Lebesgue measure on the torus T. 
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Identification of the equation. That any limiting point of the measure-valued 
sequence {Q'^'^}n>o is supported on the weak solutions of (|2.10|l follows if the 
event corresponding to the violation of the limiting equation has probability zero 
in the limit. 

For each test function <j) E C^'^([0, T] x «) that arc differentiable in time and 
twice differentiable in space, and each finite time T we define the function 



(3.8) 



1 ^ / • \ 1 ^ / ■ \ 



i=l 
N 



i=l 



t 1 



^^9s0(s, ^)x^,aAr(s)fis + i ^ ^ (/)" ( S, ^ ) a(Xi,aA' («)) 



X6 '^{Xi-bN,cN{s) - 2Xi^cN{s) + Xi+bN,cNis))ds, 



and the event 



(3.9) 0,V,= sup \V{t)\>e 

, 0<t<T 



We will prove in the next sections that for each e > we have 



(3.10) limsup P^:UOt^b,c,e)=0- 

a,h,c — — >oo 



The proof of the result (|3.10|l is complicated and is divided into several steps. 
The function that defines this event can be written as a sum of functions, see the 
beginning of section^ We will deal separately with each function in the sum and 
show that it converges to zero in probability. 



4. Identification of the limiting equation 

In this section we establish that the event 1)3. 9|) is negligible in the limit, and 
hence any weak limit Qt is supported on the solutions of the Cauchy problem 
(I2.1U|I . To save space, we are suppressing the time dependence of the test function 
(j) that defines the event H3.9I) . 

Note on notations. Assume / is some local function. We denote by Av^I^'^ jT^/ 
the average of shifts of /, namely 



21 + 1 
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We write V{t) 
(4.1) 



(4.2) 

(4.3) 
(4.4) 



N 



1 



N 



2Jo iV^1 



N 



N 



N ■ 



i - 1 
N 



N 



i + 1 



N 



— >w,ds - 



1 



1 



N 



(4.5) 
(4.6) 



i=l 
N 

E 

i=l 

2j/ ^^E'^'V^ 



7V2(w, - Av;+'/_,^i(;,)ds + 



1 



1 



2 7o ^ 



ds + 

7V2ioo(Av;+'/_;^rVr)ds 



-d{xi^aN)b '^{xi-bN,cN — 2xi^cN + Xi+bN,cN) 



ds. 



We proceed to prove that each term in the sum above converges to in probability. 

Martingale estimate (the term (4.1)). We call MN{t) the term (4.1). From 
Ito formula we know that the process {Mj^{t)}t>Q is a martingale and 



T ^ N 





•J i=i 



i - 1 
N 



- 2(j) 



N 



i + l 
N 



n 2 



aids 



The test function ip is chosen to have continuous second derivatives, therefore 
£;"<=q[M|,(r)] is of order 0{^). We can use the Doob's inequality 



pneq 



sup \MN{t)\ >e\< \e[MI{T% e > 0, 

0<t<T 



to conclude that the martingale is negligible in the limit, i.e. 



N- 



lim Pj;^'=^<^ sup |Mjv(OI > e ^ = 0, e>0 



0<t<T 



The term (4.2). A straightforward computation involving the Chebyshev in- 
equality and the entropy inequality H4.1(J|I proves that the term (4.2) converges in 
probability to 0. The test function needs to have continuous fourth derivative. 



A technical Lemma. We shall prove a Lemma that reduces the problem of 
establishing the negligibility of an event to finding that the largest eigenvalue of a 
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Schrodinger operator is negative. For an operator A : Ti ^ TL acting on a Hilbert 
space Ti. we denote by supspec^A the largest value in the spectrum of A. 

Lemma 4.1. Let {x{t)}t>o be the slope process with generator H^.^p . Under the 
assumption 



(4.7) 



limsupsupspeCi2(^-i) ( ay + — ) = 



lim sup sup 

N^oc p,E<--q[p2] = l 



aE'^^igp^ 



77 



Dn{p) 



< 0, a 7^0 



it follows that the event | g{x{s))ds > e| has negligible probability or: 



(4.8) 



lim P^"^ 

iV— *oo ' 



g{x{t))dt 



>(\=Q, e>0. 



Proof. We can use Chebyshev inequality to reduce the proof of (|4.8|) to 



(4.9) 



lim 



g{x{t))dt 



0. 



Since we do not have much information about the initial noncquilibrium distribution 
we use the entropy inequality to replace the noncquilibrium distribution in (|4.9|l by 
the equilibrium distribution. 

Before we continue, we remind that given two probability measures v and /i on 
the same probability space such that v is absolutely continuous with respect to we 



define the relative entropy of v with respect to /i by H{v\^) — ^ log ^ 



dfj. 



where 



^ is the Radon-Nikodym derivative of v relative to /i. The entropy H{iy\ij), always 
a positive quantity, is the optimal constant that makes the entropy inequality 



(4.10) 



EM] < ^|ff(^^|M)+logi?p[e"^]|. 



true for any bounded, measurable function / and a > 0. A trivial consequence of 
the entropy inequality H4.1()(l helps us to estimate i^(A), where A is some event, 

\og{2) + H{,y\fi) 



(4.11) 



iy{A)< 



In our context we use the entropy inequality for the distribution of the pro- 
cess started in nonequilibrium and the distribution of the process started in the 
equilibrium. We have. 



^ncq 
1 



g{x{t))dt 



1 



- Na ^ ^ ^ ^' Na ^ 



Na 
exp I a 



Nag{x{t))dt 



Ng{x{t))dt 



< 



< 



< - + -^logE'"i 



exp iVa / g {x (t)) dt ]+ cxp [ - Na g{x{t))dt 



a Na 

As a consequence of the inequality 

log(a + 6) < max(log(2a),log(25)) < log(2) + max(log(a), log(6)) 
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for two positive numbers a and 5, H4.9|l follows as soon as we have 



lim sup — log E°'^ 



exp Na / g{x{t))dt 



< 0, a ^ 0. 



A trivial consequence of Feynman-Kac proves our Lemma, 



(4.12) 



1 



exp Na / g{x{t))dt 



iV4 



< TsupspeCi2(^^q)(a5 + —Ln). 



N 



□ 

The microscopic current w can be replaced by local average of cur- 
rents (the term (4.3)). We shall show that the term (4.3) converges to zero in 
probability. As a consequence the current w is replaced by a local average of w, 
that is closer to a deterministic value, w by itself is a single fluctuating random 
variable. 

We check that the hypothesis of Lemma 14.11 is valid for the function gi^N — 
X^i^i {wi — Av^jt\_iWj) . Recall that if the test function has continuous 

second-order derivative the quantity 



Avf' ,0(A)_0(J_) 



N ' 



fi+l ^ 

^ N 



21 + 1 



,2 

is of order 0{4t2). Now, on integrating by parts twice, it follows for a fixed function 



/9, with E'=i[p^] = 1, that 



N 



N 



'N 



iV2 



Therefore, 



5] V^X,(p)V^(Av}tt;0(|^) - 

E^-(Av-_,^(i^) 

1=1 

{E^ip'^El^ (Av;+t,0(^)-0(^))^.]) 



< 



lim sup sup 
and hence (|4.7|) is satisfied. Moreover, 



Dn{p) < C—^—Dn{p). 



N 



Dn{p) 



iV4 



lim P'''"^ 

l,N—^oc ' 



TV 



i=i 



> e ^ = 0, e > 0. 



Inserting the fluctuations (the term (4.5)). Let /(a;_s, . . . , Xg, xq,;) G 
C^(R^*+^) be a local function that depends on the slope configuration in a box 
of size s and on the mean slope in a box of a large size I. We use the notation 
h — I — Vl for a slightly smaller I. We want to show that the fluctuations approach 
zero, in the limit or that 



(4.13) lim 



LN- 



ncq 
N,T 



N_ 

N 



N . . s 



> e ^ = 0, e > 0. 
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We apply Ito formula, 



(4.14) 



1 



N -N^ 



N , . 



(Av;tt;,^V)(^(r)) - (Av;t^L,^rV)(a:(0)) 



The summand in (|4.14|) converges to zero, as the function / is bounded. The 
second part of H4.14|l approaches zero because the L? norm of ^^j^i'^'' is of order 
O(^), as we can see below. Let 5 = ^ ^^^^ 0( J^)(Av}tti,^V)(a^), then 



E 



ncq 



21 



7V4 



_ ^ncq 



i^NQ - "^gLNgdt 



iV2 



^ncq 







_ ^ncq 


[i 







1=1 



The function kY^tk-h'^^ f 

is a cylinder function that depends just on the sites n 



such that k~l<n<k + l. Therefore the vector field Xi is zero when acting on 



most of the summands inside Av'^'j^^ , f. There are no more than 21 sites k such 



that Xi{Avjl^j}_i_^T^ f) 7^ 0. We put all these arguments together to conclude that 

Replacing the current by the Laplacian of the slope field (the term 
(4.4)). In our model the instantaneous current, w, can not be written as the 
discrete Laplacian Th — 2h — T~^h of some local function h, thus we use the method 
of Varadhan |17| for computing the hydrodynamic scaling limit of our model. The 
main idea is that the current decomposes as 

(4.15) w ^ a{xi)Ax + L^f 

for a suitable coefficient a(x;), where xi is the average slope in a cube centered at 
the origin of microscopic side I. The equation (|4.15|) is known in the literature as 
the fluctuation-dissipation equations. We have explained before that terms of the 
form Loof have no effect on the macroscopic scale. A new feature is characteristic 
to our model due to the complex interaction of the system: after filtering off the 
fluctuations from the current we are left with the Laplacian of some function and 
not a gradient, as happened for models previously considered in the literature. The 
precise meaning of (|4.15|l is given below. 

Our aim is to prove the existence of a sequence {fr}r>o of local functions and 
of the transport coefficient a such that 



limsup P^x) 



{i/:fsH^)[A 

I— 1 ^ ' 

oo/r]di| > e| 



(4.16) 



-i—li 



a{xu)Ay';JUi^^)i 



- Av''+'i 



The local function fr{x-s, ■ ■ ■ ,a;s,xo,;) depends on two arguments, the mean slope 
xo,i in a box of size I and the slope field inside a box of size s, the size s being 
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much smaller than I. We can assume that the operator ioo does not act on the 
first argument xo,i, since we can show that the action of the operator Loo at the 
boundary sites is negligible. To be more precise, Av*^'^^_;^T-'Loo/r stands for 

2/1 + 1 ■ 

Note that the function Av*^'-^_;^T-' Loo/r depends just on the value of the field inside 
the box centered at i and of size I. 

As before we use Lemma 14.11 to conclude that the event (|4.16|l has negligible 
probability in the limit. An additional difficulty shows up. If Ac is the largest 
eigenvalue of the perturbation, L + eW, of a negative operator L with principal 
eigenvalue 0, the eigenvalue Ag has the formal series expansion, 

Ae = + eE^W] + e'^ <W, {-L)-^W >^ +0{t'), 

Hence if the potential W has mean zero, one expects limg^o AfC^^ W, {—L)^^W > 
Fortunately for suitable potential W the central limit variance < W, {—L)^^W >i, 
converges to zero. 

We shall need in our context, a particular result about the largest eigenvalue of 
a perturbation operator, whose proof is found in Quastel [?]. 

Lemma 4.2. Let W be a real potential that satisfies 

< u, Wu >^< l-'^/^Di{uy/^\\u\\2, I < 

for some C small enough, or 

||M^IU<C, K<iCe)-'/^ 

Provided that the generator I'^Li has spectral gap of order one the following estimate 
holds 

(4.17) e-'^r^supspecL2(^^){l^Li - el^W) < I < W, {-Li)-^W >^ 

Spectral gap. Indeed the generator I'^Li of our model, defined by 12.6|l . has a 
spectral gap of order 1. The proof is standard by the method of Bakry-Emery (see 
Chang and Yau or Deuchel and Stroock [2]). 

The operator Li is an unbounded operator defined on the subspace C5"(R^'+^) 
of the Hilbert space L'^{v^^i) and is negative definite, with spectrum included in 
the negative semiaxis of the real line. is an eigenvalue of the operator Li but the 
eigenspace corresponding to this eigenvalue is quite large, being infinite dimensional. 

It is not hard to see that we can write the Hilbert space V^iy^^j) as the direct 
sum ®y(zK2 L'^{i'y i)- Moreover because of the ergodicity of the dynamics (|2.6() 
on the level sets of the function j/o,; = iyoiiVod^ know that the eigenspace 
corresponding to zero of the restriction of the operator L/ onto each Hilbert subspace 
LP'iVy is one dimensional. The next eigenvalue of the restriction L/|^2(t^c is a 
negative number. The distance between the largest eigenvalue and the next largest 
eigenvalue of the operator ii|L2(yc ^ is called the spectral gap of the operator, 
because it is the gap in the spectrum of the operator. 

Lemma 4.3. (Spectral gap) There is a universal constant C that does not depend 
on the conserved quantities j/o.i — {Uq i - Vq i) G such that 

(4.18) ^E,^Jp']<<{-Li)p,p>,.^^ 
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for any mean-zero function p, Eiyc Jp] =0. 

Since the generator Li of our model has a spectral gap, see Lemma 14.31 we can 
introduce the central limit theorem variance in our context. 

Definition 4.1. Suppose we have a cylinder function . . . Xs) such that 

Eijc ^ [/] =0 for all possible values of y € . Recall that Vy ^ is the canonical 
measure in a box centered at and size s defined in section\^ The central limit 
theorem variance of f on the box Ai^i is defined to be: 

(4.19) Vu(/,Z/) - 2(20 < Av;tt;,^V,(-i.,0"'(Av;tt,^rV) >.j,„, • 

// the box Aij is centered at then we use the shorter notation V; (/,?/) for the 
CLT-variance. At this point we stress that Vi_;(/, y) is a local function depending 
on the field inside of Ai^i, more precisely depending on the conserved quantities 
yi,i — {yliTyfi)! the mean slope and the linear mean of the slope field. 

At this point we stress that "Vij{f,y) is a local function depending on the field 
inside of Aij, more precisely depending on the conserved quantities yij — {yl i, yfi), 
the mean slope and the linear mean of the slope field. 

The strategy is to give a bound for the largest eigenvalue of a perturbation 
operator in terms of the CLT-variance Yi^i{f, y). Extra care must be taken because 
the CLT-variance ^i,i{f,y) is uniformly- in-y small on bounded sets and not on 
unbounded sets. 

The canonical measure ^ for our model has been obtained by condition- 
ing the grand canonical measure on the configurations with fixed mean slope and 
fixed linear mean slope in a box centered at i and of size I. The second condi- 
tioning makes the canonical measure not to have identical marginals. Actually 
the expected values of the marginals depends linearly on the site. However the 
finite-dimensional marginals of the canonical measure converges toward the finite- 
dimensional marginals of the grand canonical distribution as the size of the box 
approaches infinity, see Lemma [5m To benefit of this fact we will replace the CLT- 
variance V;(/, y) in a box of size I with its expectation E^'^\Vi{f^ y)\yk] with respect 
to the canonical measure j, in a box of larger size k. We let k go first to infinity. 
We formalize below. 

Let us define the function g as 
(4.20) 

A^^Uw, - a(xM)Av;tt JA^), - A^fXij' L^U 
Thanks to Lemma [4. II the event (|4.16|l has negligible probability if 
limsup supspec£,2(j,<!q)(g -I- 2/3— -Ljy) < 0, /3 > 0. 

r,l,N—>oo ^ -'^ 

We write the operator g + 2j3^LN as a sum of operators and we estimate the size 
of the principal eigenvalue of each operator in the sum. 

(4.21) 5-^2/3— Lw = f7i+ 1^2 + f^3 + f^4 



N 



1 
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where, 



ni = 



^2 = 



N 



^ ^ f I \^ JSl"^ 



N 



BN ^ 

i—l 



Ar4 



N 



N 



Yl'f^\J^j E'''^\V^^l{w - a{xi^i)Ax - Loo/r, ?;)|y»,fc]l|y, ,|<5. 



The operators U,2^ f^s and are understood as multiphcation operators. 

The operator Oi. Assume that M is an upper bound for the test function |(/)|. 
Define 

We have 

supspeCi2(^oq^(f]i) < 



< sup sup E'"^ 
|A|<j\/p,£;'='i[p2]=i . 



XN'^gojp'^ - ^V/(u) - a{xi)Ax - L^ofr, y)p^ 



< sup sup sup 

|A|<M!y,GR2 lp,E'=^[p^\yi] = l 

— ^V;(w - a(xi)Ax - Loo/r, y) 

P 



E'^^[\N^g,,p']^^^D,.^{p) 



Integrating by parts we can show that there is a constant C, not depending on I 
and the values of the conserved quantities, such that for each density p: 



E''^[gQ,iP^\yi]<C 



Hence the hypothesis of Lemma [4. 21 is satisfied and 



hm sup sup 

N^ao p,E'"i[p'^\yi\ = l 



E^^[\N^g,,p^]~^^D^.y) 



1 



{Zl + L)l N^oo p,E''i[p^\yi] = l 



A2 



' E^'^[^^^^^±^g,,p^]- pi'D..Jp^) 



< 



< —Yi{w - a{xi)Ax - Locfr,y)- 
P 
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The convergence on the hne above is uniform over the set of aU possible values of 
the conserved quantity yi = {yj, yf) £ M^, therefore the principal eigenvalue of the 
operator fli becomes negative as iV — > oo. 
The operator 5^2 • Let us call 

Vi.k = Vij{w - a{xij)Ax - Loofr,y) - E^'^[Vi^i{w - a{x^^i)Ax - Loofr,v)\v%,k\- 
We observe that 0,2 contains Vi^k without being multiplied with a factor N"^ . Then, 

, „ , , , ^ Ar4 

SUpspeCi2(^^q)(f22) = SUpSpeCi2(^,^q) 



N 



< 



< sup sup 



sup 

p,E-'i[p^\yk] = l 



2 



2k + 1 



sup sup N"^ 


sup 




|A|<A/afceR2 







A2(2fc + l)fc4 2 

— — ""'^'P 



Moreover, Vi^k is a bounded function (see Lemma [5.2|l and E°'^[vQ_k\yk\ = 0. We 
can apply Lemma 14.21 

"A2(2/c + l)fc4 



lim sup 


sup 


£;cq 


Af— >oo 


.p,£;<=<l[p2|yfe] = l 





-'"o,kP 



< 0. 



The convergence on the line above is uniform over all possible values of the conserved 
quantities yi — {yl,yf); therefore limsupjy^oo supspec^2(-j,^)(ri2) < 0. 

The operator fl^. We refer to Lemma l4.ll Rather than proving that 
Iimsup5 ,, ;_fc_jv^^supspeci2(^^=i)(f^3) = we will show 



(4.22) 
where 



lim sup — log E'^'^ 

S,r,Lk,N^oo N 



exp' 



ds = 0, 



O3 = 



f3N 



N 



E'^^lVii'w - a{xi)5x - Loofr,y)\yk]'^\yk\>s- 



The equation l|4.22|l follows from the estimations 



N 



logE""^ 



exp N I ri^ds 



<^logE'^^ 



1 

TJo 



exp{NTn3)ds 



N 



Uf E'^^'i^iiw - a{xi)Ax - L^fr.y)\yk]l 



\Vk\>S 



< 



< logE,.^[exp(M^E'''i[Vi{w - d{xi)Ax - L^fr, y)\yk]l\y,\>5])]- 

The dominated convergence theorem can be applied (see the note at the end of 
section O and gives us, 

lim logE,.^ [exp(M2s^q[V,(i« - d{xi)Ax - L^fr, y)\yk]l\y,\>s]) ■ 

The operator flii. The main purpose of sections [3 is to show that there exists a 
sequence of functions {fr}r>o such that for any 6 > 0. 

limsup sup E'^'^lYiiw - d{xi)Ax - Loofr,y)\yk] "^0, and 

r,l,k—>oa \yk\<S 
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sup sup E°'^[Yi{w - a{xi) Ax - Loo fr,y)\yk] < oo. 

The term (4.6) It follows that the term (4.6) converges to zero in probability 
from the following two lemmas. We refer the reader to Bertini, OUa and Landim 
PI or Guo, Papanicolau and Varadhan [7| for the proof. The proof uses mainly the 
entropy inequality H4.1()|l and consequence of Feynman-Kac formula H4.12|l . 

Lemma 4.4. (Local ergodicity) Let f be a cylinder function. Define / : R — > R to 
be the function f{a) = E^s^ [/]. Then for any 5>0, 0:T— >]Ra smooth function, 



t 

lim sup Pj 



jneq 
N,T 



N 

1 \ ^ , / I 



N^'^(nI " /(Av}tt.^,(s))) 



ds> S\ =0. 



Lemma 4.5. (Two-Block estimate) For any continuous function g : R ^ M, let 

Fk,a,N = I Av^iAv}t^^,^(5(Avj;t,2;,(.)) - 9{Ay=it--M^)))''^' ^ ^ 
Then for any 5 > 0, 

limsup PN''j-{Fk,a,N} ^ ^- 
k — *C30,a — *0, — ►CJO 

5. Computation of the central limit theorem variances 
In this section we compute the value of the limit: 

limsup sup E^g'^[¥i{f,y)] 

i— >oo |a|<(5 

for a particular class of cylinder functions / to be described later. 

We shall need a result that relates canonical and grand canonical measures, 
known as the equivalence of ensemble. The equivalence of ensemble says that 
asymptotically the marginal in a fixed box of the canonical measure is the marginal 
of the grand canonical measure. 

Lemma 5.1. (Equivalence of ensemble) Let f{x-s,...,Xs) be a bounded, local 
function. Then, for any e > 0, there exist Ni G N and Si > such that 

\E,.Jf]-E,.j-[f]\<e 

as long as k > Ni, I > Ni and — a| < Si, yf. £ R. 

We shall need a new notation, namely for the adjoint of the vector field Xq = 
d-i — 2do + di with respect to the inner-product of the Hilbert space L^{a dv'^j^). 
The adjoint is given by the formula, 

X*{h) = -Xo{ah) - (V'ix^i) - 2V\xo) + V'{xi))ah. 

Note that the current w is equal to -^^^0(2), and the slope Laplacian Aa; is equal to 

Lemma 5.2. Let h(x^s, . . . , x^) be a bounded, cylinder function such that f — 
XqIH) g L^{v^^) for all a G R. Then there exists C{h) < 00 that depends just on 
the function h such that 

(5.1) sup 1<{-U)-H E ^'/)' E ^'f>^i,<C{h). 



i,yo,i 



|j|<i-VI |'t|</-VI 
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In addition if the function h is hounded away from zero, h > C > 0, we have the 
lower bound, Ci{h) > that depends just on h, 

(5.2) Ci(/i)< sup \<{-U)-\ E ^'/)' E ^V>.s,- 

Proof.^ Call Y,{f,y) = \ < {-Li)-HEui<i-Vi^' f)^^H\<i-Vi^' f >^iy We use 
the variational formula 

(5.3) V,(/,2/) = sup— 



ID,^ , (u) 



We have 



iii<i-yi b\<i-vi 



b\<i-Vi 



1 2 



< 2lC{h)D^.^{u), 



therefore the upper bound for the variance is established. 

For the lower bound, we may chose a particular function u in namely 
u = Y^^^_l \xi. Note that Xi{u) = 1 for any —l + l<i<l—l. It follows, 

1 2 



;i?,c (u) 



7l<;-i '^iJ 



>Ci(/i). 



We have just used that a, and /i are functions bounded away from zero, and a is a 
bounded function. 

□ 

Note. Lemma 15.21 is true for any local function / = X]i=-s where the 

functions hj are bounded, local functions. 

Definition 5.1. For a bounded, local functions f such that E^c ^[f^ — 0, for all 
possible values ofyER^, we define the semi-norm: 

1 



{{f))i^ hmsup -E,. 



I y-^ 



<{-u)-\ E -V), E -v>.^, 

\i\<i-Vi |j|<'-VI 

We saw in Lemma IF?^ that {{f))a is a finite number as long as / is equal to (/i), 
where /i is a bounded, local function. By polarization we can extend the semi-norm 
(( ))q to a semi- inner product (( , ))a- For the remaining part of this section 
we compute the value of the semi- norm {{f))a for certain functions /. 

Lemma 5.3. Assume that g is a bounded cylinder function with bounded first 
derivatives, f = Lo^g, w = Xo{a) — {V'{x-i) — 2V'{xo) + V'{xi))a and Ax ~ 
X-i — 2xo + xi, then the following identities hold, 

2n 



a) {{L^g))l = 



a(a;_i,a;o,a;i)( Xo{ J2jez'^^9 



b) ^ 4:E^g'=[a{x-i,xo,xi)], 

c) {{Loo9,w))a ^ 2E^g-=[a{x-i,xo,xi)Xo{ J2jei.'^^ 9)] 

d) {{L^g,Ax))c. = 0, 
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e) {{w,Ax))a = 4. 

Proof. One checks directly using equivalence of ensemble lemma 15.11 and the 
asymptotic shift invariance of j, that the relations a)-e) hold. In particular, for 
d) and e), it is important to notice that Aa; = x^i — 2xo + xi = X^{-). 

□ 

Definition 5.2. Define the Hilhert space Tia to he the closed linear span in L^(a dv^) 
of the function 1 and functions = Xq (X^jez '''"'s) ' 'where g is a bounded local func- 
tion with bounded first derivatives. 

It is not hard to see that if / is equal to Xq(/i) then, 



{{f,Locg))a = E^gc 



iProjnA2h)Xo(j2^- 



, {{f,w))^ = E,.^[aProjnA^m- 



On both lines above Projua stands for the projection opertor in the subspace TLa- 
We are left to calculate {{f))a for a function / = XQ{h). As we will show in the 
following lemma, ((/>>« = E^s^^[a{ProjnA2h)Y]. 

Lemma 5.4. Suppose f is equal to Xq(/i), where h is a bounded cylinder function. 
The value of the semi-norm of f , 

{{f))a=E,.MProjnA2h)n 
Proof. Using Cauchy-Schwartz inequality it follows that. 



1 



((/))„= hminf -E,. 

l^oo.yt, — >a I 



> 



\i\<i~^/i \j\<i-Vi 
> E,.^^[a{PromA'^h)f]. 
Consider g = Y.\j\<i-^i'r^ f ^ where / = X^h) then 

J<(-^i) ( ^/)' ''f>-^i,= '^^ip 1 

\t\<i-Vi \3\<i-Vi p^°"I.Sp'^ 



2P 



2/2 



a{x^j,Xj,Xj)Xj{piyh 

■\3\<l~Vl 

Above pi is some maximizer function with D^i ^(fi) — 21. 
Call ui = TlT.\■J\<l-^l^'AXJPl)■ 



E,.^^[a{ulf] < E^go 



jD,g^{pi)<2. 



^ Y: {r-AX,pd? 

The above inequality shows that {ui}i is a bounded sequence in L'^{adv^) and 
hence has a weakly convergent subsequence in L?'{adv^). Let u be the weak limit 
of a convergent subsequence of {w;};. It follows that 



limsup —Ely 



I — >oo,y^ — >a 



I 



< 



\i\<i-Vi lil<'-VI 



< 
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< 2E^g'=[auh]'^ = 2E^9c[auProjH^hY < 
< 2E,g^[au^]E,g^[a{Projn^hf] < E,g^[a{Pro]H^2hf]. 
The key point that has allowed us to write the above inequalities is that the 
function u has the property Xa{T^u) ~ Xi,{t"'u) for all integers a and b, and a 
function with this property belongs to Ha, (see Lemma [2. II or Savu ^Sl)- The 
next estimate shows that Xair'^u) ~ Xb{T°'u) is valid in a weak sense. Consider a 
smooth test function 0, with bounded first derivatives. Assume a > b. We have 

<X„(r%),0>.g— < 4 E r''-\pi),X;x:cf>>,g^, 

lil<;-VI 

and 

< XaiT^Ul) - Xbir'^ui), ^ >,ge = 



a+l-Vl 



21 



a-l-i + \/I 



21 



By Cauchy-Schwartz inequality, we obtain 

<X,(TS)-^fc(T"U/),(/<>2ge< 

- / ^ a-{-l — y/l ^ a— 1 — /+VT 

- r^Xb-M) + Ji ^ r^Xb-,{p. 

j=b+l+l-Vl j=b-l + Vl 

-1-i+VI i-Vi 



< E„ 



21 

C{a - b) 



{21) 



j=h-a-i+\/I i=h-a+l+l--/l 



(5.4) 



< 



C(a - b) 



C(a - b) 



As I converges to infinity, the sequence {w;}/ approaches u in the weak sense. 
Combining this fact with (|5.4|l we can establish Xa{T^u) = Xbir^-u). 

Now we can conclude that if we have a local function / such that / = A'g(/i) 
then 



((/,/»«= lim ]e^.^ 



\i\<i-Vi \j\<i-Vi 
E,g^[a{Pro]nA'2h)f]. 



□ 

The results proved so far in this section allow us to conclude that if 5 is a local 
function equal to X^ih) and b{y\) is a coefficient that depends on the mean-slope 
in a box of size k then 

{{g + bw + Loof,g + bw + Loof))a= lim E^c ^\Vi{g + bw + Loof ,y)] ^ 

I — >oo.y^ — >a ^' 



Ea.g^ [{Prom^ {2h) + 2b + Xo{J2 ^7))']- 
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Since Proj-H^{2h) is a function in the closed linear span of the functions 1 and 
-^oiJ^iez'''^ f) where / is a local function it follows the existence of a sequence of 
functions {/r}r>o SLud of a coefficient b = i such that 

lim {{Ax + bw + L^fr))a = 0. 

r — *oo 

Lemma 5.5. Let a* > be an upper bound for the function a. For any e > and 
S > there exists a smooth function . . . ,xi,a) such that 

sup {{w — £00/ — a{a)Ax))a < e, and 

\a\<S 

sup {{w — Loaf — a{a)Ax))a < 2a* . 

a 

Proof. We calculate 

{{w + L^f - a{a)Ax))^ = "^E^.^Hl + Xo(^t7))2] - a{a) 

Let us introduce the notation Ti.b for the set of cylinder functions f{x-i, . . . , xi) 
with I I/I loo < B and ||a,/||oo < B for ^ = -/,..., L Let A{f,a) = E,.'4a{l + 
^oCigz'^V))^] - and Ai^sia) = inf/g^,^ a). The function Ai^sia) 
is upper semicontinuous and nonincreasing in I and -B, and for each a G R, 
linij^B^oo A^sia) = 0, therefore, 

lim sup Ai^sict) = 0. 

Therefore we can find for each a a cylinder function /(a) such that A{f{a), a) < e, 
if |a| < 5 + 1. We extend / to be zero on |q!| > S + 1. Then on |q!| > S + 1 we 
have that A(/(a), a) = Ej^e'=[a{x_i, xq, Xi)] — a{a) and yl(/(a), a) < a* . The avoid 
the problem that / is not smooth, we take the convolution of / with a smoothing 
kernel (p. The required function is the convolution / * 0. For a complete argument 
see Lemma 2.6 in Quastel |13| . 

□ 

Note. From the equivalence of ensemble lemma I^TI we know that for any bounded 
cylinder function / and any e > 0, there exit Ni E N and 61 > such that 

lE'^^lYiiw - aix^,i)Ax ^ L^f,y)\yk] - {{w - d{a)Ax - L^f))^\ < e, 

as long as fc > iVi, I > Ni, and |?/^ — q;| < Si. Then Lcmma l?31 helps us to conclude 
that for any e > there exists a bounded cylinder function / such that 

limsup sup E°^[Vi{w - a{xi)Ax - Locf,y)\yk] < '^f-, and 

l,k^oo \yk\<S 

sup E^'iiYiiw - a{xi)Ax - L^fr,y)\yk] < 3a* . 

a* on the line above is the upper bound for the function a in the hypothesis of 
Lemma 15.51 

Properties of the transport coefHcient. The transport coefficient a has 
been defined as the unique real number such that 

(5.5) inf ((u.-a(a)Aa;-Loo/))' =0. 

It is important to notice that the transport coefhcient is also given by the formula: 
{{w,Ax))a _ E,.42la] _ 1 



a{a) 



{{Ax,Ax))o. E,.^[{ProjnAl)ra] E,<,^[{ProjnAiWa] 
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= inf 



(5.6) 



a 



a{x-i, xo, xi)(l + Xo{2_^T^ g)y 



6. HyDRODYNAMIC limit of the model for SURFACE ELECTROMIGRATION 



In this section we will show that the model for surface electromigration |2 . 81 has 
a hydrodynamic scaling limit as well. We will prove Theorem 12.21 The proof is 
inspired by Quastel |13| . From Cameron-Martin-Girsanov formula we can find the 
Radon-Nikodym derivative 



jpneq 

"■^N,E,T 

"a/" t 



'■^^^\\^ [ ^(^' Vo'ixi-i, Xi, Xi+i)dB., 



~\Y1J^ E'^(t,J;^^o■{x^-l,Xi,Xi+l)dtJ 
and the relative entropy 



gX^y E'^yt,—)a{xi^i,x,,x,+ijdt 



Hence there exists a constant C such that H{P^'^ tI^n^t) — CN . We can use the 
inequality (|4.11|) and the superexponential estimates ^6.'A\ . (|3.4|) to conclude that 
the sequence of probability measures t\n>q is tight in the topology of the 

space X. Denote by Qe.t a weak limit of a subsequence of {Q^\ t\n>o- 

We proceed further to identify the limit Qe,t- To save space we ignore the time 
dependence of the test function cj). 



^ Ell 4>{j,)xm - ^ Ef=i 4>{j,)xm - 



T ^ N 



2N ■ 



— J2 ( N'w, + E(t,-)a,) r{^)dt + M^iT) 



(6.1) 
(6.2) 



2N ~ a{xi^k){^x)i - Loot'' Jr)(t>" { — )dt + Mn{T) 

i=l 



2N 



N 

E 

N 



E[t,-] (a, - a{x,^k)) + N^{LooT'fr - LN,ET'fr) 



^"{j^)dt 



N\^x), + E[t,- 



<t>"{j^)dt 



Above {MAr(t)}t>o is a martingale and as in section 4 we can prove that is negligible 
in the limit. Because of the entropy inequality H4.1()|l and the superexponential 
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estimate 



lim sup — log E*^^ 

r,LN^oo N 



exp 



<o, 



the event on line 1)4.16(1 is negligible under P^]ij rp. For the new model that we 
have the negligible fluctuations are L^^eI and not Loo/ (i-e., the term (6.1) has 
no contribution towards the limit of the model). The contribution coming from 
nontrivial fluctuations Lm^eJ — L^of are gathered in the coefficient in front of 
the vector field E in the nonlinear equation (|2.14|) . Because of this reason, the 
coefficient turns out to be the transport coefficient d, defined by the variational 
formula (|2.11|l and not a{a) = E^g'^ [a] 

We shall show that the sequence of smooth local functions {/r}r>o introduced in 
the Note right after the Lemma 15.51 can be chosen to have the additional property 
that the term (6.2) is negligible in the limit. Note that the term (6.2) is negligible 
if limsup^ fc -P^s.tI 1^(^)1 > e} = where, 



(6.3) 

We can write 
(6.4) F{T) = 

(6.5) 

(6.6) 



I 

N 



^ 1 ^ f i \ 



N 



N 



N 



I 1 ( i\ 



>s 



Recall that in the proof of Lemma [5.51 we introduced A{f,a) to be the difference 
E^,s<= [a(l + -^o(X]iGZ ■''V))^] ~ a{a). We will show that we can modify the sequence 
{fr}r>o to have the properties 

(6.7) sup A(/r(a), a) < a*, sup lE^a^ [ago] — a{a)\ < 5a*, 

Q,r a,r 



(6.8) 



lim sup A(/r(Q!), a) = 0, and 

'■^°°|a|<5 



(6.9) 



lim sup |£'j,9<: [ago] ~ a{a)\ = 0. 



We write, 

Eug" [ago] - a(a) = {E^g- [ago] - E^g^^ [ag^]) + {E^gc [ag^] - a(a)) 



= E^g<= [ago{go - 1)] + A{fr{a), a). 
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Since |i?i/g<= [ago(5o ^ 5: 4a* and A{fr{a), a) < a* uniform in a and r, (|6.7(l 
follows. Consider f*{a) to be the minimizer of A{f,a) among local functions of 
x-i through xi. f* has the property that for any fr, 



E„ 



0. 



Now, thanks to Cauchy-Schwartz inequality, 



lE^g-iagoigo - 1)]| = 



< 



< E,.^[a{Xo{J2r^fr)n^'E,.^[a{XoiJ2T^ir - fr)))V. 
We modify fr such that 



sup E'^gc 

\a\<S 



a(Xo{J2r^r~fr)) 



< 



6a* 



To conclude (6.7) is negligible because of local ergodicity lemma lOI (6.8) is negli- 
gible because for any fixed S we have ()6.9|l and (6.9) shrinks to zero once S becomes 
very large. Also it is important to have property (|6.7|l . Our result follows. 

□ 
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